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ABSTRACT 


This paper presents a new numerical method for solving 
Schroedinger's equation in radial form in order to get a set 
of wave functions for use as starting values in the solution 
of the Fock equation. As a necessary adjunct to the solution, 
an improved numerical method is given for calculating potentials 
for use in the Schroedinger equation. Finally, a method is 
outlined for solving the Fock equation numerically, using the 


wave functions obtained above as starting values. 


fe 
jute 





ENERGY OF THE MOST LOOSELY BOUND SLECTRON 


FOR IONS OF COPPER, SILVER, AND GOLD 


If the energy and distance from the nucleus of the most 
loosely bound electron in an ion of a metallic lattice can be 
found, it is possible to use this data to calculate certain 
physical constants of the metal, e.g. the elastic constants. 

Wigner and Seitz! have shown how to find this binding 
energy, neglecting exchange forces, by solving Schroedinger's 


equation in radial form: 


gq, -# &R .[vir)-ele = 0 


2m ay 


for various values of E. It is found that for each such value 


of E there are several radii at which the condition 


AR 
(2) ar 


Kk 
C 


is satisfied. Most of these points are very close to the 
nucleus and are essentially the same for all values of E, If 
the outer two radii are plotted against E a curve results 
which exhibits a minimum, It is this minimum which gives the 
required data, 

Fuchs“ has extended the procedure to take into account 
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exchange forces. His method involves solving the Fock equation: 


'p, Wigner and Seitz, F., Phys. Rev., 43, 804 (1933) 
eK, Muchs, Proc. Roy.ssocasmA, 15158505 195m) 


28, Fock, Zeits., F. Physik, 81, 195 (1933) 





ee) 


(3) — a+ [vi)- ER - \ lor elrJAe ae Clee 


The most practical method for solving this equation is 
first to solve equation (1) and then use the wave functions 
obtained to evaluate the integral in (4). This makes the right 
side a function of the independent variable instead of zero. 
When this resulting equation is solved the new wave functions 
are then used to re-evaluate the integral and the procedure is 
repeated as required. 

As the first approximation, then, a solution of (1) is 
needed. To get it, it is necessary first to calculate the 
values of V to be used. 

From the work of Hartree and ners tables of values of 
Z, the effective charge, are available as functions of C » the 
radius measured from the surface of the nucleus. These values 
have been accepted in this paper, although there may be errors 
in them since most were calculated by hand before the era of 
the electronic computer. The required values of v, however, 


were recalculated, where 
(4) | 


and 2, is the effective nuclear charge for potential which is 


different from Z. The values of v were calculated by 


4p, Hartree, Proc. Roy. Soc., A, 141, 282 (1933) 


A, Douglas, Hartree D., and Runciman, W. Proc. Cambridge 
Phil, Soc., 51, 486 (1955) 
M. Black, Manchester Memoirs, 79, 29 (1934) 
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integrating the field intensity 
a a ( 
f= Zr) uh 
— t pee -_ ee Ro r 
(5) ) ak (2) a w( ye oe , 


The usual method? of calculating v is to solve the 


differential equation 





Zep tpo-F 
(5) —— = —f 
2G % 


by Euler's method, which is simple but inaccurate, It can be 


ee that the error is of the order of 


e 


Some simple algebra and substitution reduces this expression to 
(3) e: LOG) \2(@) -z(e,)| 
aaa 


Since Z is assumed to be linear between any two adjacent radii. 
Craph 1 gives the relative % error for Cut resulting from this 
calculation. It is large enough so that a more sophisticated 
method of solution, equation (5), was decided upon. The 
evaluation was done using a Gaussian 3-point quadrature formula, 
which is based on the assumption that between any two neighboring 
radii Zis approximated by a third degree polynomial. 


To evaluate the integral in (5), 
Cs. f a 
(9) A( ¢,) = \ [ ui 


aD, Hartree, Proc. Cambridge Phil. Soc., 24, 114 (1928) 


oxunz, K., Numerical Analysis, McGraw-Hill, N.Y., 1957,pgs 177-8 





it is first necessary to normalize the limits to the range O to 1 
in order to use readily available values of the Gaussian coef= 


ficients, This can be done by the substitution 


, f20-(4+0) | 
(10) w= FZ ( . oe - 


Equation (9) then becomes 


> ( (2u- \(e@- P.) + (exe) ) 


an ale): 46-2 eanecer eeer 


The coefficients for the Gauss J=-point formula are 


5 
ae Wi, = iP 
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Thus, (11) becomes 


15) Ale) Aé2|5§ 


where 


oo) 2 (Fer%es) 


— 
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(KC, +0 &)” ( CPi +e C.)” 








\-3 z (Gx ) 


(¢,+ Cr)* 


(14) «<= \.17454¢666 5 Qe 222540 334 
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It can be showm' that the error is of the order of 


(15) 625 — I$ (e+ S(ed|-Al@)- ae Wits $url 


where f represents the integrand in (9). The magnitude of the 





relative error for Cu’ is plotted in Graph 1. It is never 


greater than .25% 


To use (13) it is necessary to develop an interpolation 


‘Lanezos, C., Applied Analysis, Prentice Hall, N.J., 1956, 
pes 404-7 
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formula for calculating the 3 numerators on the right from the 
available table of values of Z. This has been done using a 
Gregory-Newton second degree divided difference interpolation 
formula, A third degree formula which would be consonant with 
the basic assumption involved in the Gaussian 3=-point method 
would have increased the labor of programming significantly 
and would have made no difference in the final result, to the 
limits of accuracy retained. 


The resulting three por ae are of the form 


2(e)- 2(()+ = eee [ele.)-2(¢,)] 


(15) ee é\(< C+ee, Fal 
eset: (o)-elew-Leraeler, 


Appendix 1 gives the CDC~1604 program for calculating v, 





starting at = =6 Bohr radii and proceeding inward. v( é) is 
calculated by equation (42) for C =5.5 Bohr radii as 
~.18181818181818, The first v calculated by the program is 
that for C. =>) Bohr radii. The values obtained are given in 
Table 1, 

Equation (1) can now be solved. It is first desireable to 
write it in terms of Bohr radii and Rydberg energy units to cet 
a dimensionless equation, 


Let 
Va 


(17) C= a, 


where a. is a Bohr radius. Then 


(18) aA*R 1 Atk 
ae Gee (C= 
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so that (1) becomes 
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where € is a dimensionless number, since a Rydberg unit is 














given by 
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where 


(27) Q=2-2r 


Equation (26) can be solved readily using a Runge-Kutta 


4th order approximation. Let 
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After making the calculations (32) and inserting in 


formulas (41) we ger 
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Cot, : : 
As before, Q (RE calculated using a Gregory- 
Newton second degree divided difference interpolation formula 
of the same form as equation (16). 

Since equations (34) give only successive increments in 
Rand y it is necessary to obtain starting values for each. 

Assume that near the origin 
(35) Re | bot b, ‘i ~be~ +. 
Also, for small @ , 
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_ eC 
where Z is the complete charge on the ion. Differentiating 
(55) twice and inserting it and ™ in (26) gives 
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Inserting (39) in (35) 


so that b, = |. 
As starting values, let us calculate initial values for 


Rand y at © =.002. Inserting the coefficients in (5) 
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We can ignore © compared with 22°, Then the starting 


values for any € are 
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The CDC-1504 program for the calculation is in Appendix 2. 
The magnitude of the initial energy was set at .4 in each case 
and increased in increments of .01 to 1.0 (to 1.6 for gold). 
The potentials for @ =.002 were calculated from equation (36). 


Potentials for C>6 were calculated from 
\ 
(42) 7 


The wave functions from this calculation are intended 
merely for starting values in equation (3), so no great accuracy 


: : . 5 
1s required, The error® is of the order of h . Only a sample 


Sxunz, Ke, op. cit., pe 188 
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curve of KR vs fe is shown for each metal, that for the energy 
value at which the minimum of the curve of € vs C was achieved, 
This is given in Craph 2. The curves of € vs e are shown in 
Cran 9. 


To solve (4) it may first be Fon in dimensionless forn: 
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Here the fy 1 refer to the wave functions for the individual 
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electronic shells. These are available (cf note 4). 
For Cut, the filled shells are the 1s, 2s, 2p, 3s, 5p, and 


oad, corresponding to n,l values of 10, 20, 21, 30, 5i55%2.ne ie 


integral then becomes (using r< @ as an example) 
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and similarly for r > 
All the integrals on the right can be evaluated using 
numerical methods with the wave functions calculated above as 


starting values, although the programming will be tedious 
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The interration, for all practical purposes, probably need not 
proceed beyond the range 0 to 6. The right hand side of (46) 


then becomes 


(47) Fi) =a F,deleat (arash (ear (A 
is 55 (p)+ a +, (e) 


where the a's are the constants resulting from the integrations. 


It is then necessary to solve the set of equations 


Te | &(¢)- é|R+F(e) 
ie* tt 


by a Runge-Kutta approximation using equations (31) and (32) 


ke 
hi 


(48) 


as before. The resulting wave functions are then used to 
re-evaluate the intecrals on the right of (46) and so on. The 
procedure is terminated when succeeding sets of wave functions 
differ by less than the required degree of accuracy. 
Unfortunately, there was not enough time available to write 


the CDC-1604 program. 


Ut 





CONCLUSIOINS 


The results obtained differ sienificantly from thoseof other 
investicators. For example, Fuchs calculated the minimum energy for 
copper as approximately .8, occurring at about 2.9 Bohr radii. 

He gives no estinate of his possible error, however. 

Since the error in finding each increment in the method 
used here is of the order of h, the cumulative error in R to the 
radius at which the minimum occurs is in no case greater than ,003. 
The accuracy of the method cannot be called into question, therefore. 

There is always the possibility of a mistale in algebra or 
programming, sucn as the omission of a factor of tuo. Every step 
in the calculation has been checked over four times, and if there is 
such an error it must be very subtle in nature. 

It should not be thought, however, that the results of 
previous investifators in this field are necessarily correct. They 
were calculating by hand, and the methods which they were forced to 
use were so oversimplified that reasonable accuracy was often 
impossible. Indeed, if the previous investigators (who all seen 
to have used the same methods) had been in error only .01 units 


for each increment of R, their total error could be very large. 
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APPENDICES 


6: Fileures in the 2° eolurn refer to Aas 
aie S8 i + : 
and in the 34 column to Au’, where different 


COlumm, olen as rer Ct 





ANS 


VR1 


SFT 


ORG 
EQU 
EQU 
EQU 


LDA 


LDA 
BMU 
OTA 
LDA 
FMU 
FAD 
OTA 
FDV 
STA 
LDA 
FLU 
ofA 
LDA 
FLU 
OLA 


LDA 


APPENDIX I 


10000 


11000 


12000 


11100 


VIN 


VR1 


00000 


00000 


ALPHA 


RAD+1 


STOR 


BETA 


RAD+2 


oTOR 


STOR 


TWO 


STOR+1 


STOR 


STOR 


STOR 


BETA 


RAD+1 


STOR+2 


ALPHA 


RAD+2 
20 


. - 


SALAS | 


“re 
poaedey wei 





FAD 0 STOR+2 
STA O STOR+2 
FDV O TWO 
STA O STOP+3 
LDA O STOR+2 
FMU 0 STOR+2 
STA O STOR+2 
LDA 1 RAD+1 
PAD 1 RAD+2 
STA O STOR+4 
FDY O TWO 
STA 0 STOR+5 
LDA O STOR+4 
FU =O STOR+4 
STA O STOR+4 
LDA 1 ZETA 
FSB 1 ZETA+1 
STA OO STOR+S 
LDA. 1 ZETA+1 
FSB 1 ZETA+2 
STA 0 STOR+7 
LDA 1 RAD 
FSB 1 RAD+1 
STA 0 STOR+8D 
LDA 1 RAD+1 
FSB 1 RAD+2 
STA O STOR+9D 





LDA 


FDV 
PMU 
FAD 
FAD 
FDV 


STA 


So CFO. Cleo owe =>. ©. Oo Oo 2S eS) 


oO OF 90 Oo oO -—- 


— 


O 


RAD 
RAD+2 
STOR+10D 
STOR+6 
STOR+9D 
oTOR+11D 
STOR+/ 
STOR+8D 
STOR+12D 
STOR+11D 
STOR+12D 
STOR+11D 
STOR+ED 
STOR+9D 
STOR+1 OD 
STOR+9D 
STOR+1 
RAD 
STOR+15D 
STOR+8D 
STOR+5 
STOR+1 4D 
STOR+1 
P*D+1 
STOR+13D 
STOR+10D 
STOR+11D 
STOR+14D 
ZETA 
pLOR 


STOR+1 


faye 





FDV 


FLU 


OLA 


LDA 


FSB 


FDV 


EMU 


FAD 


FAD 


FDV 


PAD 


FU 


STOR+4 
RAD 
STOR+15D 
STOR+8D 
STOR+5 
STOR+16D 
STOR+3 
RAD+1 
STOR+15D 
STOR+10D 
STOR+11D 
STOR+16D 
ZETA 
STOR+2 
STOR+1 


FIVE 


STOR+3 
STOR+5 
RAD 
SLOR+1752 
STOR+8D 
STOR+5 
STOR+18D 
SleR+5 


RAD+1 








i reper al 







CONV 


FMU 


FDV 


FLU 


FAD 


FDV 


PMU 


FAD 


PMU 


FDV 


STA 


LDA 


FMU 


FAD 


SLI 


0 


* 


y 


STOR+17D 
STOR+10D 
STOR+11D 
STOR+18D 
ZETA 
STOR+4 
FIGHT 
STOR+3 
TWO 

NINE 
STOR+3 
RAD+2 
RAD+1 
STOR+3 
VR1 

VR1 
RAD+2 
ANS 
ANS+1 


00002 


00000 


4 ANS 


h 


{1000 


ISK 


ah 


1 


29D 


ISK 


1 


Lop 





E‘I 0 00000 
001 O ANS 
00 
INI 4 00006 
4. 00000 
SAU O CONV+1 
ENI © 00000 
ISK 3 17D ISK 4% 19D ISK> 3) 16D 
SLJ 0 CONV 
ENI © 00000 
SL3 0 00000 
VIN DEC -.18181818181818 
ALPHA DEC 1.774596666 
BETA DEC 22540334 
THO DEC ay 
FIVE DEC 5 
BIGHT DEC oe 
NINE DEC 9 


RAD DEC Ge 


DEC 


2 

D 
DEC 45 

h 

3 


DEC : 

DEC 8 
DEC 3.6 
DEC 54 


2) 





DEC 


DEC 


DEC 


DEC 


DEC 


BC 


pee 
De 

260 
2,6 
2.4 


ae 


1.8 
1.6 
1.4 
1.3 
ee 


ee 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


1.4 
ee 
1.1 


ue) 


2 ® 
Oo \O 


An ~ 


DEC 


DEC 


DEC 


DEC 


UPS. 
1.4 
Le, 
1.2 


ibe 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DiC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


Dac 


DEC 


DEC 


DiC 


16 
14 


212 


0025 
woe 


015 


S000 
.006 


0044 


ay 


ZETA 


DEG 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


uC 
| 
p=) 


DEC 


DEC 


DEC 





DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DisC 


15050 
1.068 
1.095 
1541350 
1.179 
1,243 
1.550 
1.448 
1.606 
Lee? 
O07 
2.507 
SOO 
5.896 
4455 
5.148 
6.018 
Teite 
8.457 
10.119 
12.051 
14.153 
15.196 
16.193 
17.114 
17.964 


ZETA 


DEC 


DiC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


Dec 


DEC 


ee 


14,21 
15.94 
These 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DiC 


DEC 


DEC 


V2 
aye 
1.95 
Pero 
aoe 
eet 
444 
D625 
Oe 
(Bue 
8.88 
10.18 
11.68 
13.34 
1519 
Liat 
192s 
21513 
23.34 
25.032 
27.82 
30.96 
34,73 
38.94 
10.65 
42.36 





18.791 
19.736 
20.187 
20.699 
21.279 
21.929 
22,646 
2515 
24.209 
SAG 
25.689 
26.520 
26 .508 
26.884 
om 
27.451 
Bia eo 
28.106 
28.280 
28,451 
28.615 
28.759 
20 5C15 
26.948 
28.992 


29.000 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


END 


Ab O5 
45,62 
47.13 
48,62 
50.09 
Digee 
53.40 
55.44 
Dito 


Orv 
Oo 

e 
\N 
a \ 


62.92 
65.41 
67.76 
70.14 
72.86 
15093 
11.29 
19. 





ae Pa 


sia LPL ae 
rit A 


eee io0ny 
‘ise 


. 





4 





DEG 45.5 
DEC 45.78 
DEC 46.06 


DEC 416,30 


DEC 46.63 
Dec 46.85 
DEC 46,91 
DEC 46.95 


DEC 46.99 


pec 47. 
Dace. 
END 


Note 1: Library subroutine DECOF 4 must be placed in memory 
starting at cell 71000. Output will be on magnetic 


tape unit 4, 


Mote 2: Program prints ¢ and VU. Two to a line in decimal 


Llcatvane pom. 





BEGIN 


RETUR 


Sik 


LDA 


STA 


LDA 


APPaNDIX II 


SCHR EQ BY RK METHOD 
20000 
22000 
21000 


21600 


INEM 


apn 
ENERGY 


>» OOOOO 


OOOO00 


OO000 


OOO0O00 


RIN 
a a 


ARL 


Ye 


WYE 


RAD 


RAD+1 


STOR 


RAD 


RAD+2 


STOR+1 


RAD 1 


RADA2 





FDV 


FDY 


STA 


LDA 


RAD 
RAD 
STCR+2 
POT +1 
POT +2 
STOR+2 
FOUR 
STOR+1 
STOR+34 
Eos 
POT+1 
oT OR+2. 
STOR+2 
FOUR 
STOR 
STOR+1 
STOR+4 
POT 
STOR+4 
STOR+3 
STOR+5 
RAD+1 
RAD 
STOR+5 
STOR+5 


ST OR+ / 





PMU 


LDA 


FSE 


ora 


LDA 


FLU 


FLU 


PRU 


FDV 


LDA 


FLU 


LDA 


PLU 


al 
on 


PiU 


STA 


LDA 


STOR+6 
STOR+85D 
POT 
ENERGY 
STOR+9D 
STOR+5 
ENERGY 
STOR+°D 
STOR+8D 
ARE 
FOUR 


STOR+10D 


ENERGY 


Wy P 


tJ 


STOR+/ 
STOR+11D 
STOR+5 
WYE 
STOR+/ 
STOR+1 2D 
STOR+5 
ARE 
STOR+S 
TWO 
STOR+14D 


ENSRGY 


\N 
WN 





Pictu 


Fil 


OTA 


LDA 


FLU 


FAD 


LDA 


FLU 


FDV 


THREES 
ARE 
STOR+5 
STOR+14D 
POT 

ARE 
STOR+5 
STOR+15D 
WYE 

SIX 
STOR+15D 
STOR+14D 
STOR+13D 
STOR+ 12D 
STOR+11D 
STOR+10D 
STOR+6 
SIX 
STOR+16D 
POT +1 
WYE 
STOR+5 
STOR+18D 
WYE 


STOR+5 





FLU 


FU 


STA 


LDA 


ARS 
STCR+19D 
STOR+5 
BNERGY 
STOR+20D 
POT +1 
EMERGY 
STOR+20D 
STOR+1 9D 
STOR+/ 
TWO 


STOR+21D 


STOR+22D 
POT +1 
ARE 
STOR+24D 
POT 

Bu BRGY 
STOR+24D 
STOR+5 
ENERGY 
STOR+24D 


ARE 





rhiU 


BDV 


Sok 


PLU 


BLU 


PiU 


FAD 


AD 


FAD 


HAD 


STOR+/ 
THO 


STOR+25D 


STOR+5 
ARE 
FOUR 


STOR+26D 


STOR+22D 
STOR+26D 
3TOR+27D 
STOR+26D 
STOR+25D 
STOR+25D 
STOR+21D 


3TO2R4+18D 


el 
© 





CORY 


FLU 


FDV 


FAD 


ols 


ep 
-] 
ad 


LDA 
STA 
LDA 


FAD 


ENMI 
001 
004 
ini 


BTA 


2 VALUE+9 
1 RAD+1 

2 VALUE 
O STOR+1S5D 
O ARE 

O ARg 

2 VALUE+1 
1 RAD+1 

2 VALUE+2 
2 00004 

O OQ0000 

1 56D 


RETURN 


_- oOo 
© 
© 
© 
O 
O 


4. 71000 


00 
5 00004 
5 00000 


QO CONV +1 


ISK 





ENI © 00000 
Isk 4 56D Isk 4 59D ISK 4 52D 
So ame Ol 


LDA O BNERGY 


bowe 9 00079 ISK 4 00170 


Tien 42pse «6-4 


DO DEC 2. 


Six DEC 6, 

DsLEN DEC -,.O1 

RIN [B)O] Oe (81872 

Yi DEC 887 DiC | DEC ./09 


RAD DsC 002 


DEG? ~.005 DEC .004 mae .02 
DEC .01 18) KOM) 0.6) DEC .04 
DEC .015 DEC .008 BEC .06 
BECe 402 Dae 0 DEC .08 
DEC .025 DEG 015 DEC) 21 

DEC .0% Dace. 02 Wee ae 





DEC 


DEC 


DEC 


214 


DEC 


DEC 


oti 


. 144 


Wes) 
XO 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


DEC 


~ 144 


1a 
1.2 


1 





ANS 


~29000 r 


FOr 


DEC 


DiC 


1.4 
eS, 
1.6 
Tyo 


ome 
2.4 
2.6 


one 


Dae 
44 


3.8 
ae 
45 
De 
Deo 





DEC 


DiC 


DEC 


=11336. 


5539 4 
a0 lon 
2648.8 


2074.6 


1694, 1 


1423.9 


1222 


945 2 

799.16 
629.14 
Bene 

458 ,28 
399.42 
312.54 
25am 

208.14 
175.07 
149.49 


DEC 


DEG 


DEC 


DEC 


DEC 


DEC 


4 


to 
27820. 
-14880. 
- 9974, 
7410. 
- 5838. 
- 5890. 
- 2876, 
~ 2258. 
- 1841, 


- 1545, 


Ds 


Duc 


DEC 


DEC 


DEC 


DIC 


DEC 


DiC 


DEC 





9.01 


6.902 
De29 


DIC 


DEC 





DEC 


DEC 


Note 1: Note 1 to Appendix I also applies here 


mote 2: Program prints ?, RE ama AR 


along each line. 


Pas 


o 


Le 


in decimal floating point 
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most loosely bound electro 


wu KAT 


3 2768 0 
DUDLEY ie ——- yy 





